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Abstract. An optimal control problem for semilinear parabolic partial differential equations is 
considered. The control variable appears in the leading term of the equation. Necessary conditions 
for optimal controls are established by the method of homogenizing spike variation. Results for 
problems with state constraints are also stated. 
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1. Introduction. We will give necessary conditions of optimal controls for parabolic partial 
differential equation (PDE, for short) with leading term containing controls. This is an analogue 
of the result we got for elliptic PDE with controls in the leading term ([11]). Let us consider the 
following controlled parabolic PDE of divergence form: 

dtz(t, x) — V- (A(t, x, u(t, x))Vz(t, x)j = f(t, x, z(t, x), u(t, x)), in £1t, 
z(t,x)=0, on[0,T]x5fi, (1.1) 
z(0, x) = zq(x), in fl, 

where SIt = (0, T) x £1, T > and £1 C R" is a bounded domain with a smooth boundary 
dfl, A : J1t x U — > R nx ™ is a map taking values in the set of all positive definite matrices, 
/ : x R x U R, with U being a separable metric space and zq(-) £ -L°°(f2). Function u(-), 
called a control, is taken from the set 

U = {w : fir — > U | w(-) is measurable }. 

Under some mild conditions, corresponding to a «(•) £ IA, (1.1) admits a unique weak solution 
z(-) = z( - ;u(-)) which is called the state. We measure the performance of the control by the 
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following cost functional 

J(u(-))= f°(t,x,z(t,x),u(t,x))dtdx (1.2) 
Jn T 

for some given map /° : Ot x R x [/ -> E. Our optimal control problem is stated as follows. 
Problem (C). Find a u(-) e U such that 

J(«(-))= inf /(«(■))■ (1-3) 

Any u(-) e satisfying (1.3) is called an optimal control, and the corresponding f (•) = z(- ; 
is called an optimal state. The pair (z(-), u(-)) is called an optimal pair. When A(t, x, u) = A(t, x), 
Problem (C) has been studied by many authors, see [10] and the references cited therein. Works 
concerning the elliptic cases with leading term containing controls can be founded in [4] , [5] , [6] , 
[7], [11], [12], [13], and [14], etc. However, it seems that there are only few works devoted to 
parabolic cases (see [3], [16], etc.). 

In this paper, we make the following assumptions. 

(51) Let T > and f2 be a bounded domain in R" with a smooth boundary dfl. 

(52) Let U be a separable metric space. 

(53) Functions A(t, x, v) = (ay (t, x, v)) take values in the set <S™ of n x n (symmetric) positive 
definite matrices, which are measurable in (t, x) G VLt and continuous in v e U. Moreover, there 
exist A > A > such that for almost all (t, x) G fir, 

A|£| 2 <(A(t,x, «)£,£) <A|£| 2 , VCeR", veil, (1.4) 

where ( • , • ) stands for the inner product in R™ . 

(54) Functions f(t,x,z,v) is measurable in (i, x), continuous in (z,v) € R x U, and continu- 
ously differentiable in z. Moreover, there exits a constant M > such that 

zf(t,x,z,v) < M(z 2 + 1), V(t, x, z, v) € £It xExf (1.5) 

and for any R > 0, there exists an Mr > such that 

\f(t, x, z,v)\ + \f z (t, x, z,v)\ < Mr, a.e. (t,x,v) e O t x [/, |z| < i?. (1.6) 

(55) Function f°(t,x,z,v) is measurable in (t,x), continuous in (z,w) € R x U, and continu- 
ously differentiable in z. Moreover, for any R > 0, there exists a Kr > such that 

\f{t,x,z,v)\ + \f°(t,x,z,u)\ <K R , &.e. (t,x,v)e{l T xU,\z\<R. (1.7) 



Our main result is the following. 
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Theorem 1.1. Let (S1)-(S5) hold and z € L°°(Q). Let (z(-),u(-)) be an optimal pair of 
Problem (C). Let ^S(-) be the weak solution of the following adjoint equation 

d t tp(t, x) + V- (A(t, x, u(t, x))Vtp(t, x)) = f°(t, x, z(t, x),u(t, x)) 

-f z (t,x,z(t,x),u(t,x))ip{t,x), in fl T , 
•ip(t,x)=0, on[0,T}xdn, 



tp{T,x)=0, in fi. 



The 



H(t, x, z(t, x),ip(t, x), S7z(t 7 x),Vtp(t, x),u(t, x)) — H(t, x, z(t, x), tp(t, x),Vz(t, x), Vi[)(t, x), v) 
> ^\A(t,x, v)~^(A(t,x,u(t,x)) - A(t,x,v))Vz(t,x)\ \A(t,x,v)~^ (A(t,x,u(t,x)) - A(t,x,v))Vip(t,x)\ 

+ i ( A(t,x 7 vyi(A(t,x, u(t,x)) - A(t,x, v))Vz(t,x), A(t,x, vyi(A(t,x,u(t,x)) - A(t,x,v))Vi>(t,x)), 

Vw e U, a.e. (t,x) € fir, (1.9) 

where 

H(t, x, z, ip, r},v) = { tp, f(t, x, z, v) ) -f°(t, x, z,v)-{ A(t, x, v)£, n ) , 

(t, x, z, ^(,i),»)e[0,T]x!lxRxEx R" xR"x[f. (1.10) 

Since the right hand side of (1.9) is always nonnegative, (1.9) implies 

H(t, x, z(t, x), %j)(t, x),Vz(t, x), VV>(t, x), u(t, x)) 
= maxHtt, x, z(t, x), tp(t, x), Vz(t, x), Vr/i(t, x), v) , 

y v e u, a.e. (t,x) e n T . (1.11) 

When A(t,x,v) = A(t,x), the right hand side of (1.9) is zero, thus, the result automatically 
recovers those for the classical semilinear case without state constraints ([10]). 

Since U is not necessarily convex, it is well-known that people usually use spike variations to 
derive necessary conditions for optimal controls. Such a spike variation technique does not di- 
rectly work for problems with leading term containing the control. To overcome the difficulty, we 
adopt the idea of homogenization for PDEs to carefully select some special type spike variations 
of controls so that we can have desired "differentiability" of the state with respect to the control. 
We can see in [11] that such a method is useful for the cases of elliptic PDEs. The main idea to 
treat parabolic case is same to that for elliptic case. However, there are some new difficulties in 
studying properties of variational equations. 

Comparing Theorem 1.1 and the corresponding result for elliptic case in [11], we can see that 
they are similar when n > 2 and slightly different when n — 1. More precisely, Theorem 1.1 of 
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this paper is very similar to Theorem 1.1 in [11] for high dimensional cases. In particular, for 
parabolic case with n = 1, instead of 

H(t, x, z(t, x), tp(t, x),z x (t, x),i> x (t, x),u(t, x)) - H(t, x, z(t, x), ip(t, x),z x (t, x),$ x (t, x),v) 

(A(t, x, u(t, x)) — A(t, x, v)) 2 . T . , , . . — /, 

> V V K .'' / " z x (t,x)^ x (t,x), VveU, a.e. (t,x) e Sl T , (1.12) 

./II X, U ) 

we have (1.9), i.e., 

H(t, x, z(t, x), tp(t, x),z x (t, x),ip x (t, x),u(t, x)) - H(t, x, z(t, x), ip(t, x),z x (t, x),ip x (t, x),v) 
(A(t, x, u(t, x)) - A(t, x, v)f 



> 



, ~ x (t,x)ip x (t,x) , V v e U, a.e. (t, x) e fl T - (1-13) 

One can see that (1.12) is similar to the corresponding result for elliptic case with n = 1, while 
(1.9) (i.e., (1.13)) is similar to the corresponding result for elliptic case with n > 2. We mention 
that for elliptic cases with n > 2, the corresponding right hand of (1.9) follows from a fact given 
in Lemma 2.5. While for parabolic case with n = 1, the right hand of (1.9) (i.e., (1.13)) follows 
in a different way. In fact, it follows from (1.12) and 

H (i, x, z(t, x), ip(t, x),z x (t, x),$ x (t, x),u(t, x)) - H (t, x, z(t, x), ip(t, x),z x (t, x),$ x {t, x),v) 
> 0, Vuef/, a.e. (t,x) G T - (1-14) 

From the proof of Theorem 1.1, one can see that (1.12) can be yielded from using spike variation 
along space-direction and (1.14) can be yielded from using spike variation along time-direction 
(see (3.24)). 

Another difference between parabolic cases and elliptic cases appear in that there are three 
possible types of homogenized equations for parabolic cases when taking a different scale for 
the time and the space variables, while there is only one type of homogenized equations for 
elliptic cases. Difficulty occurs in analyzing the second type of homogenized equations (see the 
proof of Lemma 2.2 for details). Despite the different types of homogenized equations, the 
variational equations are same and we finally get same optimality conditions for the three cases. 
Nevertheless, we think results of this paper will be useful to analyze the second-order variational 
equations, which is a problem more difficult than that for first-order variational equations. 

The rest of the paper is organized as follows. In Section 2, we present some preliminary 
results. Section 3 is devoted to a proof of our main result. Problem with state constraints will 
be discussed in Section 4. 

2. Preliminaries. In this section, we will give some preliminary results needed in proving 
Theorem 1.1. For Y = [0, a\] x [0, a 2 ] x . . . [0, a n ], a function g(x) on R" is called F-periodic if 
it admits period aj in the direction Xj (j = 1, 2, . . . , n). 
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Lemma 2.1. Let r > 0, S E (0, 1) and (SI) hold. Let h(-) E L 2 (n T ), and A m {-) = (a™(-)) e 
Z,°°(Q T ;<S!J) sucA that for some A > A > 0, 



A|£| 2 <(.4 m (M)£,0<A|£| 5 



V£ E R", (i, x) e T , m = 1, 2, 3, 4. 



(2.1) 



G(t,x,s,y) = (gij{t,x,s,y)) = (g lj (t,x,s,y 1 )) = < 



Define 

' A\t,x), if ({*}, { yi })e [5, l)x [5,1), 
A 2 (t,x), if (W,{i/i})e[o,J)x 

A 3 (t,x), if ({*},{!,!})€ [5, l)x [0,5), 

A 4 (t,x), if (W,{l/i})e[0,*)x [0,5), 

(2.2) 

w/iere {a} denote the decimal part of a real number a. For e > 0, let z £ (-) E L 2 (0, T; Hq(Q)) be 
the weak solution of 



d t z £ (t, x) - V- [G(t, x, 1 ^) Vz £ (t, x)] = h(t, x), in T , 



(2.3) 



(2.4) 



(2.5) 



z £ (i,x)=0, on[0,T]xffi, 
z £ (0, x) = z (x), in il 

with z (-) E L°°((l). Then 

z E {-)^z{-) 1 weakly in £ 2 (0, T; Hq(Q)) 
with z(-) being the weak solution of 

dtz(t, x) — V- (Q(t, x)Vz(t, x)) = h(t, x), in Or, 
z(i,a;)=0, on [0,T] x 90, 
z(0, x) = 20(2)5 in 
and Q(-) = (qij(-)) E L°°(Qt',S+) being given by 

Qij{t,x)= / dyi / (gij(t,x,s,yi) + gn(t,x, s,yi)d yi ip> (t,x, s,y!))ds , l<i,j<n, (2.6) 
Jo Jo 

where ^ fe (i,x, •) e L|(0, 1; T4^' 2 (0, 1)/R), # means the function is [0, 1] periodic. 
For r < 2, cp k (t,x, ■) is the unique solution of 

d Vl (gik{t,x,s,yi) + gii(t,x,s,y 1 )dy 1 tp k (t,x,s,yi)) = 0. 
For r — 2, <p k (t,x 7 •) is the solution of 

d s f k (t, x, s, yi) - d Vl (gik{t, x, s,yi) + gn(t, x, s, yi)d Vl ip k {t, x, s, y x )) = 0. 
While for r > 2, f k (t,x,s,y) = if k (t,x,yi) is the solution of 

dyi (Jo 9lk ( t,X,S,Vl " >ds + / 9ij( t , x ,s,yi)dsd yi tp k (t,x,yi)^) = 0. 



(2.7) 



(2.8) 



(2.9) 
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Proof. The above proposition is a corollary of Theorem 2.1 in [2], Chapter 2 (see also 
Remark 1.1 and "Comments and Problems" there). The result can also be got by the technique 
of two scale convergence ([1], [8]). According to [2], 



Qij{t,x) 



r f 1 n 

/ d JJ {9ij(t,x,s,y) + ^)2 i g lk (t,x,s,y)dy k (p 3 {t,x 1 s 1 y))ds , l<i,j< 
J[o,i] n Jo k=1 



For r < 2, <p k (t,x, •) G i|(0, 1; VF^' 2 ((0, l) n )/R) is the unique solution of 

n 

^2 d yi{9ij{t,x,s,y)8 jk + g i:j (t,x,s,y)d V:j tp k (t,x,s,y)) = 0, 

where 5^ equals to 1 if i = j and if i 7^ j . 
For r = 2, ip k {t, x, •) is the solution of 

n 

d s (p k {t,x,s,y) - d yi(9ij{t,x,s,y)5 jk +gij{t,x,s,y)d y:l ip k (t,x,s,y)) = 0. 

While for r > 2, </5 fe (t, x, s, y) = <^(t, x, y) is the solution of 

y~] dyii / 9ij{t,x,s,y)8 :jk ds + I gij(t,x, s,y) ds d yj tp k (t,x,y)) = 0. 



(2.10) 



(2.11) 



(2.12) 



(2.13) 



Since G(i, x, s, y) is independent of y 2 , 2/3, . . . , j/ n , we must have (f k (t } x, s, y) — f k (t 7 x, s, y\) 
and consequently, (2.10)— (2.13) becomes (2.6)— (2.9). □ 



The following lemma concerns the "derivative" of Qij in 6 = 0. 

Lemma 2.2. Le£ r > 0, A > A > 0. Assume A > a r „ > A, \b m \ < A, \c m \ < A (m = 
1,2,3,4). Let S e (0, 1) and 



(a 5 ( S ,y),6 5 ( S ,2/),c 5 ( S ,y))= ^ 



(ai,6i,ci), if ({*},{!/})€ [5, l)x [5,1), 

(02,62,02), if (W.M)e [0,5) x [5,1), 

(a 3 ,6 3 ,c 3 ), if ({>},{*/})€ [5,1) x [0,5), 

(a 4 ,6 4 ,c 4 ), if ({*},{!/}) €[0,$)x [0,5). 



(2.14) 



(a3.63.c3) 


(oi,6i,ci) 


(0.4. 64, C4) 


(a2.62.c2) 
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Let </>{(■) e L|(0, 1; W^ 2 (0, 1)/R) be the solution of 



d y (b s ( S ,y)+a s (s,y)d y( l > l( S ,y)) = 0, 



^2(") be the solution of 



d s 4(s,y) - d y (b 5 (s,y) + a 5 ( S} y)d y 4( S ,y)) = 



and 4>l(s,y) = 4>l(y) be the solution of 



dy(J b 5 (s,y)ds + J a s (s,y)dsd v 4(y)) = 0. 



Then there exists a constant C = C(A, A), independent of 5 € (0, 1), such that 



^ 5< ^« A,Si w , ( fo 3 - &i)(c 3 - ci) 
c°{s,y)d v (l>i{s,y) ds + — 



5 ./(I ./() «3 

Proof. . I. It follows from (2.15) that 

6 5 («, y) + <*'(*, J/)M(*>W) =!>'(*)■ 
Since ^f(-) is [0, l] 2 periodic, we get that 

b s (s,y) 



<CVS, k = l,2, 



i s {s,y) 

Solve p s (-) from the above we get 



j' p^s) b\s,y) dy= f M{[s y)dy = ^ s£[0jl] _ 
Jo a \ s tV) Jo 



d v (t>i{s,y) = < 



(1 - 5)a 3 + Sai 
_ b 4 -b 2 5 

(1 — 5W + 5a2 
63 - 61 



(1 - 5)a 3 + 5ai 
&4 — 62 



5, (s,t/)e(5,l)x(5,l), 
(*,!/)€ (0,5) x (J, 1), 
(1-5), (s,t/)e(5,l)x(0,5), 



Thus 



(1 - 5)a 4 + 5a 2 
2A 



(*,l/)| < 4 



A direct calculation shows that 



A 



2A 



(1-5), (a,!,) e (0,5) x (0,5). 
5, ( S ,t/)e[0,l]x(5,l), 



— (1-5), («,!,)€ [0,1] x (0,5). 



+ 



(63 - 6i)(c 3 - ci) 



«3 



5(2 - 5) ( ^ 3 ~y (C3 ~ Cl) - 5(1 - 5) ^ " 52)(C4 " C2) 



(1 — 5)03 + 5ai (1 — 5)a4 + 5a 2 

4A 2 . . .. 4A 2 12A 2 . 
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II. It follows from (2.16) and the periodicity of ^O) that 

= / [\b 5 ( S ,y)d v 4( S ,y)+a 5 ( S ,y)\d y 4( S ,y)\ 2 }d S dy 
Jo Jo 

= 11 [(b 5 (s,y)-b 1 )d v 4( S ,y)+a 5 ( S , y )\d v 4( S ,y)\ 2 }dsdy. 
Jo Jo 



Thus, 



l r i 

9 v 4>i{s,y)\ 2 dsdy 
o Jo 

1 ^ ^ 



< - I I a°(s,y)\d y <f> 6 2 (s,y)\' dsdy 



o Jo 



= -\J o 1 j\b 5 ( S , y )-b 1 )d y 4( S ,y)d S dy 

\{J J \ bS ^ y) - 5 i I 2 ds d y) * ( J f \ d v^ y)\ 2 ds d y) 

i 

\d y (t> S 2{s,y)\ 2 dsdyY . 



to Jo ' v Jo Jo 

2V2AV5 ■ fl rl 



o Jo 



Therefore, 



II A ^11 <r 2 ^ A ^ 

Il c 'l/<P2|U 2 ([0,1] 2 ) < 



A 

On the other hand, (2.20) implies 

H 5 y0lllL 2 ([o,l] 2 ) ^ Y~ 
Denote $ 5 (-) = <^(-) - we have 

ii^iu, n ,,,,< 2(1 + f AV? . 

For any p(-) € W£ 2 (0,1), it follows from (2.15)— (2.16) that 



/ [b 5 ( SlV )+a 5 ( SlV )d y cj>[{ S ,y)]d yV (y)dy = ^ V*e[0,l], 
Jo 

L 

[b s (s,y) + a d (s,y)d y (j) s 2 (s,y)}dyip(y) dsdy = 0. 



l r i 

$i 



o Jo 
Thus, 



f f a S (s,y)d y $ 5 (s,y)d yl p(y)dsdy = 0. 
Jo Jo 

Therefore / a s (s, y)d y Q s (s, y) ds is a constant and consequently, 
Jo 



a 5 (s,y)d y $ d (s,y)ds 



o 

rl pi 



na 5 (s, w)d y $ 5 (s, x) ds dx 



Optimality Conditions for Parabolic PDEs 



/ / (a s (s, x) - ai)d y $ 5 (s,x) dsdx 
Jo Jo 

~ (/ / \ a5 ( s ^)- a ^ dsdx Y \\dy* s \\L>aom 

< V2°A^. 2(1 + ^ )A ^= 2(2 + ^ )A %, Vye[0,l]. (2.30) 



Nothing that (2.15) implies 

V.S/„ „,\o /<5/„ „ ,\ | „<5/„ „,\q .(5/ „,\o <<5 



[6 d (s,y)^( S ,y) + a d ( S ,y)9^(s,t/)^(s, 2 /)]rfy = 0, Vs e [0,1], (2.31) 
we get from (2.22) that 

/ / a s (s,y)d y $ 5 (s,y)d y <f> s 2 (s,y)dsdy = 0. 
Jo Jo 

Thus denote 

= b 3 - h = &4 - ^2 

(1 — 6)0,3 + Sai ' ' (1 — 5)d4 + <5&2 ' 

we have 

Wdy&WvaoAW = I / a 5 (s,2/)|^( S ,y)| 2 rf S ^ 
Jo Jo 

= - f f a 5 ( S ,y)d v $ s (s,y)d y <pi(s,y)dsdy 
Jo Jo 

= - f f a s (s,y)d y ^ S (s,y)e 1 . s 8d S dy~ f [ a s (s,y)d y <t> S ( S ,y)e 2 , s 8 dsdy 
Js Js Js Jo 

+ [ [ a 5 (s,y)d y $ 5 (s,y)6 1 , s (l-6)dsdy 
Jo Jo 

+ [ [ a s (s,y)d y <f> 5 (s,y)(6 2 , S -6 hS )(l-8)dsdy 
Jo Jo 

< f [\ 5 {s,y)\d y &{s,y)\ ■ 2 ±5dsdy + f f 
Js Jo A Jo Jo 



s 2A 

a 5 ( S ,y)d y <S> 5 ( S ,y)d S .—(l-5)dy 



f/ / a S (s,y)\d y $ 5 (s,y)\ ■ —(1 - 6) ds dy 



o Jo 



^ -fS\\d y ^\\ LW;a5) -VA+ 2{2+ f )A2 5 ■ *± s 



2A 
\ 

4A 
+ T 



VA<5( J" J" a s ( S ,y)\d y ^ s ( S ,y)\ 2 dsdyy 
6AVA X|| „ . 4(2 + v^)A» 



Thus 



\\d y $ 5 \\ LH[0 , m < -^\\d y * 5 \\l W;aS) < ^-S. (2.32) 
Then, (2.30) can be improved as 

J a 6 (s,y)d y 3> 5 (s,y)ds\ <V2AVS-^^6= ^J^ SVS, Vye[0,l]. (2.33) 
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Moreover, 



6 



< 



d y $ A (s,y)ds dy 
\(f o J o ( a iXis,i](y) + a3X(o,S)(y))d y $ S (s,y)ds dy) 2 



1 

< - 
" A 



+ 



a s (s,y)d y $ s (s,y)ds 



dy 



AWo Jo (^ a2_ai ' )x [ <5 ' 1 ]^^ + ^ a4_a3 ^( ' ,5 )^) 9y$ ' 5( - S ' y ' )ds dy ) 2 



< 



< 



A 

12h. 2 Vh 

a 2 Va 

12AVA 

a 2 7a 



Therefore, 



< 



< 



1 

6 
1 

+ 

12A 2 



i /.l 



( 4A2<5 / \d y ^{s,y)\ 2 ds)dy 

a aVa a 2 Va 

(6 3 - h)(c 3 - Ci) 



Jo 
1 



c (s,y)d y (j) 2 (s,y)dsdy + 
1 



03 



• / *c ^ ^ w , ( & 3 -&i)(c 3 - ci) 
/ c (s,y)aj / <?!> 1 (s,y)ds+ — 



a 3 



\J q J q c s ( S ,y)d y <S> S ( S ,y)d S dy 

\J J ( c5 ( s '^) _ c iX[8,i]{y) ~ C3X(o,6)(y))dy<f> S (s,y)dsdy 
(ciX[8,i\ (y) + C3X(o,8)(y))d y $ S (s, y) ds dy 



1 

6 

12A 2 
A 



6 + 



1 r l 



JO 

6 + 



• 1 r $ 



Jo 



((c 2 - ci)x[SA](y) + (c 4 - c 3 )x(o,,5)(2/))i9j / 4>' 5 (s,i/) dsdy 



x Jo 



d y $ s (s,y) ds 



< 



< 



< 



12A 2 r . 2Ay/S... _ 4|| 



A 5+ * 

12A 2 r 2AV5 8A-y/A r A 28A 2 VA 
■ 6 + 



1 A 2 

d y $ 5 (s,y) ds 







dy 



A (5 
56A 3 VA 



a 2 Va 

III. By (2.17), 



-Vs. 



-■S+'i- 

aVa <J a 2 Va 



/ b s (s,y)ds + / a t '(s,y)dsdy( l 
Jo Jo 



l(y) 



is a constant. Then similar to (2.19), we can get from the periodicity of 4>l(-) that 



d y <t>l(y) = { 



(b 3 -h)(l- 6) + (b 4 ~b 2 )6 
6(1 - <J)oi + 6 2 a 2 + (1 - 6) 2 a 3 + 6(1 - 6)a 4 ' 

(b 3 -h)(l- 6) + (b 4 -b 2 )6 

6(1 - 5)a 1 + 6 2 a 2 + (1 - 6) 2 a 3 + 6(1 - 6)a 4 



ye (6,1), 
(1-6), ye (0,6). 
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Thus 



1 ,-\ 



c s (s,y)d y (f) s 3 (s,y) ds dy 



(1-5) 


(1 - S)(bs - h) - 


h5(b 4 -b 2 ) 




(1- 


<5)(c3 - ci) + 5(c4 - c 2 ) 




5(1 - 5)m - 


h<5 2 a 2 + (l- 


-8 


2 a 3 


■f 5(1 - 5)a 4 



One can verify that 

1 r 1 * 



s Jo Jo 



/ 12A 2 . 



c S (s 1 y)d y (t>i(s 1 y) ds dy + 
(b 3 - bi)(c 3 - ci) 



(6 3 - 6i)(c 3 - ci) 



A3 



(l-J)2(6 3 _6 1 )( C3 _ Cl ) 



a 3 



5(1 - 5) ai + 5 2 a 2 + (1 - 5) 2 a 3 + 5(1 - 5)a 4 



12A 2 r 8A 3 c 20A 3 c 



(2.37) 



Combining (2.21), (2.35) and (2.37), we get (2.18) with C 



56A 3 VA 



□ 



The following result is concerned with the well-posedness and regularity of state equation 
(1.1). 

Lemma 2.3. Let (S1)-(S4) hold and z £ L°°(Q). Then for any u(-) £ U, (1.1) admits 
a unique weak solution z(-) £ L 2 (0, T; Hq (Q)) n L 00 (fly). Furthermore, there exist a constant 
K > 0, independent of u(-) £ U , such that 



II Z (-)IIl 2 (0,T;H1(^)) + |N(-)IU°°(O t ) < 

Moreover, there exists an a £ (0, 1), such that for any Q CC fir, it holds that 

\\z(-)\\ Ca{Qo) < C(Q ) 

for some constant C(Q )- 



(2.38) 



(2.39) 



Proof. The result is quite standard. We give a sketch of the proof. Fix u(-) £ U. Let 
m > 0, define 

{/(t, x, z, u), \z\ < m, 
f(t,x,-m,u), z<-m, 
f(t, x, m, u), z > m. 

For fixed z(-) £ L 2 (0 T ), let z m (-) be the solution of 



d t z m (t, x) - V- (A(t, x, u(t, x))Vz m (t, x)) = f m (t, x, z(t, x),u(t, x)), 
z m (t,x)=0, on [0,T]x9fi, 
z m (0,x) = z (x), in Q,, 



in fli 
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Then there exist a constant C m > such that 

\\z m {-)\\mv,T:Hi m + lk m (-)llL^ T) < c m . (2.40) 

Moreover, there exists an /? = fi m £ (0, 1), such that for any Q CC £!t, it holds that 

\\z m (-)\\c^ Qo ) <C m (Qo) (2.41) 

for some constant C m (Qo). 

Using (2.40) — (2.41), we can see that the map z(-) >->• z m (-) is continuous and compact from 
some ball of L 2 (Qt) to itself. Thus, Schauder fixed point theorem implies that the map has a 
fixed point Z m {-). We have Z m {-) e £ 2 (0, T; Hq(Q)) and 

d t Z m {t, x) - V- x, u(t, x))VZ rn {t, x)) = f m {t, x, Z m {t, x), u(t, x)), in ft T , 
< Z" l (t,x)=0, on[0,T]x<9O, 
Z" l (0,x) = z Q (x), in SI 

Noting that (S4) holds, we can modify the proof of Theorem 7.1 of Ch. 3 in [9] to get that 

\\Z m (.)\\ L ~ {QT) < C 

with C being independent of m. Let m > C, we see that (1.1) admits a unique weak solution 
z(-) e L 2 (0,T;H^(Q)) n L°°(n T ) and (2.38) holds. Finally, by (1.6), 

\f(t,x,z(t,x),u(t,x))\ < M K . (2.42) 

Thus, (2.39) follows from Theorem 10.1 of Ch. 3 in [9]. □ 



Lemma 2.4. Let 5 e (0, 1), r > 0, 



Cm{t, x ) 



5 ml , if ({*},{a;i})e[5,l)x[<5,l), 

5 m 2, ]£{{8},{ Xl })e[0,5)x[5,l), 

6 m3 , if ({s},{ Xl })€[S,l)x [0,5), 

^ S m4 , if ({s},{^i}) g [0,5) x [0,5), 

Then Cm(4^, — ) (m = 1,2,3,4) converges weakly to \i m in L 2 (VLt) with 



(t,x) e fi T - 



(2.43) 



m = (l-5) 2 , ^2=^ = 6(1-5), m = 5 2 . 

Proof. Such results are quite well-known and can be proved by modifying the proof of 
Riemann's Lemma. One can verify easily that for any rectangle F C fir, 



lim 

E-S-0 + 



/ XF(t,x)( 1 (-!-,—)dtdx = n 1 
Jqt £ £ Jn T 



XF(t,x) dtdx. 
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Since the set of all linear combinations of characteristic functions Xf(-) is dense in L 2 (SIt), we 
get that Ci(-r, — ) converges weakly to fix in L 2 (Sl T ). The remains are similar. 



□ 



Lemma 2.5. Let n > 2. Let £,77 G R" &e two nonzero vectors. Then 

sup £ V, = MH+^, 



|x|=l 



where denotes the transpose of a matrix E. 

The proof of above lemma is easy. See [11], for example. 



(2.44) 



3. Proof of the Main Theorem. In this section, we present a proof of our main theorem. 
The proof is divided into several steps. Let u(-) G U be an optimal control and z(-) be the 
corresponding optimal state. Let r > 0, u 2 (-), Us(-), u 4 (-) eWbe fixed. 

I. Homogenizing spike variation of the control. Let <5 e (0, 1) and e > 0. For any 

(t, x) — (t, xi,x 2 , • • • , x n ) € Sl Tl define 

u(t,x), if ({±},{*±}) G [5,l)x[5,l), 

u 2 (t,x), if({±},{*±})e[0,S)x[6,l), 

u 3 (t,x), if ({±},{*±})e[6,l)x[0,8), 

^ ^(t.s), if ({±},{*±}) G [0,5)x[0,6). 

Then u <5,£ (-) e Let z s - £ (-) be the state corresponding to u S ' £ (-), i.e., 

<9 t /< e (t, x) - V- (A(t, x, u s ' £ (t, x))Vz S ' £ (t, x)) = f(t, x, z 5 - £ (t, x), u s > £ (t, x)), 
z s ' £ (t,x)=0, on[0,T]xdSl, 



u s > £ (t,x) = < 



(3.1) 



in fix 



z 5 ' £ {0,x) =z (x), 



SI. 



By Lemma 2.3, there exists constants K > and a e (0, 1), independent of 5, e, such that 



II^HIIlW;^)) + ll* d ' £ (-)IU~(n T) < * 



and 



^' £ ||c«(Qo) < C(Q ) 



(3.2) 
(3.3) 
(3.4) 



for any Q CC with some constant C(Qq). 

By (3.3), for fixed 6 e (0,1), we can extract a subsequence (still denoted by itself) such 
that z s ' £ (-) converges to a function z s (-) weakly in L 2 (0, T; Hq (SI)) as £ ^ + . By (3.4) and 
Arzela-Ascoli's theorem, z S ' £ (-) converges uniformly to z 5 (-) in C(Qo) for any Q CC Six- Then, 
it follows easily from 

||* 5>e (-)IU~(n T ) < K 
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that z S ' £ (-) converges strongly to z s (-) in L 2 (Qt) and almost everywhere in fix- 
By (1.6) and (3.3), 

f(t,x,z s ' £ {t,x),u s ' £ {t,x)) - f(t,x,z s (t,x),u S ' £ (t,x)) < M K \z s - £ (t,x) - z s (x)\. (3.5) 

On the other hand, by Lemma 2.4, for any h € L 2 (CIt), when e — > + , 

/ /(i, x, a;),u 5 ' £ (i, ar))/i(t, x) dtdx 
Jn T 

-+ (1-5) 2 f(t,x,z s (t,x),u(t,x))h(t,x)dtdx + S(l-S) f(t,x,z s (t,x),U2(t,x))h(t,x)dtdx 
Jn T Jn T 

+5(1 — 5) / f(t,x, z s (t,x),u 3 (t,x))h(t 1 x) dtdx + S 2 f(t,x,z s (t,x),u i (t 1 x))h(t 1 x)dtdx. 
Jn T Jn T 

Combing the above with (3.5), we get that along a subsequence e —¥ + , 

f(t,x,z s ' £ (t,x),u s ' £ (t,x)) 
-> (1 - <5) 2 /(*, x, z 4 (t, x), u(t, x)) + 5(1 - 5)f(t, x, z s (t, x),u 2 (t, x)) 

+6(1 - 5)f(t,x,z s (t,x),u 3 (t,x)) + S 2 f(t,x,z s (t,x),Ui(t,x)), weakly in L 2 (VL T ). 

Let 

h s - £ (t, x) = f(t, x, z s ' £ (t, x),u s ' £ (t, x)) - (1 - 5) 2 f(t, x, z s (t, x), u(t, x)) 

-8(1 - 5)f(t, x, z 5 (t, x),u 2 (t, x)) - 5(1 - 5)f(t, x, z s (t, x),u 3 (t, x))h(t, x) 
-5 2 f(t,x,z 5 (t,x),u 4 (t,x)) 

and z S ' £ (-) be the solution of 

d t z s - £ (t, x) - V- (A(t, x, u 5 ' £ (t, x))Vz S ' £ (t, x)) = h 5 > £ (t, x), in T , 
z S ' £ (t,x) =0, on [0,T]xaO, (3.6) 
z s ' £ (0,x) = 0, in ft. 



Then 



\z 5 ' £ (T,x)\ 2 dx + X \Vz 5 ' £ (t,x)\ 2 dtdx < ~z s > £ (t,x)h s > £ (t 1 x)dtdx. (3.7) 
in Jn T Jn T 

As z 5 ' £ (-) converges strongly in L 2 (Qt), z S ' £ (-) converges strongly in L 2 (Qt) too. Consequently, 

it follows from (3.7) that 

z s - £ (-)^Q, strongly in L 2 (0,T;i#(fi)). (3.8) 

By Lemma 2.1, 

z S ' £ (-)-z 5 ' £ (-)^z s (-), weakly in L 2 (0, T; ^(ft)) (3.9) 
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with z s (-) being the weak solution of 

d t z s {t, x) - V- (Q s (t, x)Vz 5 (t, x)) = (1 - 5) 2 f(t, x, z s (t, x), u(t, x)) 

+5(1 - 5)f(t,x,z 5 (t,x),u 2 (t,x)) 
+6(1- 5)f(t,x,z s (t,x),u 3 (t,x)) 



(3.10) 

+5 2 f(t, x, z s (t, x),Ui(t, x)), in Ot, 



z s (t,x)=0, on[0,T]x<9O, 
z S (0, x) = zq(x), in Q, 

where Q s (-) = (gf,(-)) G £°°(0 T ;5!f) is given by 



(t,x) = d yi {gij{t,x,s,y{) + gf 1 (t,x,s,y 1 )d yi ip j (t,x,s,y 1 ))ds, l<i,j<n (3.11) 
Jo Jo 



{g S ij{t,x,s,yi)) = < 



(3.12) 



and 

aij(t,x,u(t,x)), if ({s},{ yi }) G [5,1) x [6,1), 

a ij {t,x,u 2 (t,x)), if ({s},{yi}) G [0,5) x [<5, 1), 

a i: ,-(i,ar,U3(i,a;)), if ({s},{ yi }) G [<S, 1) x [0,(5), 

aij(t,x,u 4 (t,x)), if ({s},{yi}) G [0,<5) x [0,6). 

For r < 2, ^ fc (t,x, •) G i|(0, 1; W#' 2 (0, 1)/R) is the unique solution of 

d yi (5ife(*i aii yi ) + #fi(i, x, s, yi )d yi ip k (t, x, s, y x )) = 0. (3.13) 
For r = 2, tp k (t,x, ■) is the solution of 

<9 s ^ fc (i,x,s, yi) - d yi (g s lk (t,x,s, yi ) + gl 1 (t,x,s,y 1 )d y ^ k (t,x,s,yi)) = 0. (3.14) 
While for r > 2, f k (t, x, s, y) — (fi k (t, x, y\) is the solution of 

d yi{J 9 S ik( t , x , s ,yi) ds + J gt {t,x,s,y 1 )dsdy 1 ip k {t,x,y 1 )^j =0. (3.15) 
Combining (3.9) with (3.8), along a subsequence, we obtain 

z S ' £ (-) -+ z s (-), weakly in L 2 (0,T;H% (Ct)). (3.16) 

Note that for fixed 6 G (0, 1), since any subsequence of z S ' £ (-) has a further subsequence converg- 
ing to the same z 5 (-) weakly in L 2 (0,T; H^(fl)), z S ' £ (-) itself must converge to z s (-) weakly in 
L 2 (0,T;H^)). 

In addition, by the optimality of u(-), we have 



J(u(-)) < lim J(u 6 ' £ (-)) = J 6 = ((l-5) 2 f°(t,x,z 5 (t,x),u(t,x)) 

+5(1 - 6)f(t, x, z 5 (t, x),u 2 (t, x)) + 5(1 - 5)f(t, x, z 5 (t, x),u 3 (t, x)) 
+5 2 f(t,x,z s (t,x),Ui(t,x))^ dtdx. (3.17) 
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II. Linearized state equation. We now would like to let 5 — » + . Denote 

z 5 (x) ~ z(x) 



Z\x) 



Then it follows from (3.10) that 



d t Z 5 (t, x) - V- (A(t, x, u(t, x))VZ 5 (t, x)) = V- (- 



(t,x) € ^T- 



\t,x) - A(t,x,u(t,x)) s 



Z s (t,x) = 0, 
Z 5 (0,x) = 0, 



Wz\t,x)) 

+ (1 - 5f I f z (t, x, z(t, x) + s(z 5 (t, x) - z(t, x)), u(t, x)) dsZ s (t, x) 
Jo 

+ (! - s ) (/(*> x , ^(*> a), u 2(t, a;)) - f(t, x, z(t, x), u(t, x))j 
+ (! - s ) (/(*. a>, ^(*> a;), "30, a;)) - /(i, ar, z(t, x), u(t, x)fj 
+S (f(t, x, z s (t, x),Ui{t, x)) - f(t, x, z(t, x), u(t, x))j , 

on [0,T] x 90, 

in O. 



in Hi 



By (S3), 



|ajj(t, a;, v)\ < A, V (i, x) e Qt; v € U; 1 < i, j < n. 
Thus, it follows from (3.11) and Lemmas 2.2 that as 5 — > + , 



1iS, x ) - aij(t,x,u(t,x)) 



converges in L°°(Q) to 



Oij(t,x) = a,ij(t,x,U2(t,x)) +aij(t,x,u 3 (t,x)) - 2aij(t,x,u(t,x)) 

[a u (t,x,u 3 (t,x)) - au(t,x,u(t,x))] [aij(t,x,u 3 (t,x)) - aij(t,x,u(t,x))] 



au(t,x,u 3 (t,x)) 



1 < i, j < n. 



(3.18) 



(3.19) 



(3.20) 



On the other hand, z <5 (-) is bounded uniformly in L 2 (0, T; Hq (O)). Thus, we can prove step- 
by-step that as S — > + , 2' (5 (-) is bounded uniformly in L 2 (0, T; iJp (0)), converges to z(-) 
strongly in L 2 (0, T; #o(0)), and Z 5 (-) converges to Z(-) weakly in L 2 (0, T; Hq(£1)) with Z(-) 
being the weak solution of 

d t Z(t, x) - V- (A(t, x, u(t, x))VZ(t, x)) = V- (Q(t, x)Vz(t, x)) 
+f z (t, x, f (t, x),u(t, x)) Z(t, x) 
+f(t, x, z(t, x), u 2 (t, x)) + f(t, x, z(t, x), u 3 (t, x)) 
-2f(t,x,z(t,x),u(t,x)), in Cl T , 

Z(t,x) = 0, on[0,T]x<9O, 

Z(0,x) = 0, in ft, 



(3.21) 
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where 



Q(t,x) = (9ij(t,x)) = A(t,x,u 2 (t,x)) + A(t,x,u 3 (t,x)) - 2A(t,x,u(t,x)) 

[A(t, x, u 3 (t, x)) - A(t, x, u(t, x))] e x e[ [A(t, x, u 3 (t, x)) - A(t, x, u(t, x))] 



(3.22) 



ej A(t,x,u 3 (t,x))ei 



From (3.17), we have 



< lim 



J 5 - J(u(-)) 



f°(t,x,z(t,x),u(t,x)) Z(t,x) 



o Jn 7 
+f°(t, x, z(t, x),u 2 {t, x)) + f°(t, x, z(t, x),u 3 (t, x)) 



-2/°(i, x, z(t, x),u(t, x)j*J dtdx. 



(3.23) 



III. Duality. Let ip(-) be the solution of the adjoint equation (1.8). Then (3.23) becomes 







f(t, x, z(t, x), u 2 (t, x)) + f°(t, x, z(t, x), u 3 (t, xj) - 2/°(i, x, z(t, x),u(t, x)) 



+ (d t ip{t, x) + V- (A(t, x, u(t, x))Vip(t, x)) + f z (t, x, z(t, x), u(t, x)) ip(t, x)) Z(t, x) 
= / f°(t, x, z(t, x), u 2 (t, x)) + f°(t, x, z(t, x), u 3 (t, x)) - 2f°(t, x, z(t, x),u(t, x)) 

+ ( - d t Z(t, x) + V- (A(t, x, u(t, x))VZ(t, x)) + f z (t, x, z(t, x),u(t, x)) Z(t, x)) $(t, x) 

- J 



dt dx 



dtdx 



f°(t, x, z(t, x), u 2 (t, x)) + f°(t, x, z(t, x), u 3 (t, x)) - 2/°(t, x, z(t, x),u(t, x)) 

-(f(t,x,z(t,x),u 2 (t,x)) + f(t,x,z(t,x),u 3 (t,x)) - 2f(t,x,z(t,x),u(t,x))j tp(t,x) 
+ {e(t,x)Vz(t,x),Vtp(t,x)) dtdx 

/ [ 2 ( /(t 

, x, z(t, x), u(t, x))tp(x) - f°(t, x, z(t, x), u(t, x)) - ( A(t, x, u(t, x))Vz(t, x), Vip(t, x) ) ) 

-(^f(t,x,z(t,x),u 2 (t,x))ip(x) - f°(t,x,z(t,x),u 2 (t,x)) - (A(t,x,U2(t,x))Wz(t,x),Wi)(t,x))^ 

- (/(*, x, z(t, x),u 3 (t, x))ip(x) - f°(t, x, z(t, x), u 3 (t, x)) - ( A(t, x, u 3 (t, x))Wz(t, x), Vtp(t, x) ) ) 
Vz(t, i) T (A(t,x,u 3 (t,x)) — A(t, x, u(t, x)))eiej (A(t, x, u 3 (t, x)) — A(t,x,u(t,x)))V'tp(t,x) 



ej A(t,x,u 3 (t,x))e! 
2H(t, x, z(t, x), r/S(t, x), Vz(t, x), Vf/>(i, x), u(t, x)) 
—H(t, x, z(t, x), ^(t, x), Vz(i, x), VV>(i, x), u 2 {t, xj) 
-H(t, x, z(t, x),ip(t, x), S7z(t, x),Vip(t, x),u 3 (t, x)) 

— &(A(t, x, u 3 (t, x)) — A(t, x, u(t, x)), A(t, x, u 3 (t, x j), Vz(t, x), Vr/>(t, x), ei) dt dx, 
where H is defined by (1.10), 

&(A,B,£,ri,v) = ^-^r|~ ' {A,B,£,ri,v) e S n x 5™ x R" x R n x R" (3.25) 



dtdx 



(3.24) 
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and S n is the set of all n x n real symmetric matrices. 

IV. Maximum condition. By a standard argument ([10]), it follows from (3.24) that 

2H(t, x, z(t, x), ip(t, x),Vz(t, x), V^(t, x), u(t, x)) 

> H{t 1 x 1 z(t, x), ^>{t, x),Vz(t, x), VV>(t, x), w) 
+H(t, x, z(t, x),ip(t, x), Vz(t, x),Vip(t, x) 1 v) 

+$(A(t, x, v) - A(t, x, u(t, x)), A(t, x, v),Vz(t, x), V?/>(i, x), ei) , 

\/v,w€U, a.e. (t,x) e Sir- (3.26) 

Further, it is easy to see that (3.26) is equivalent to the following two inequalities: 

H (t, x, z(t, x), V>(t, x), Vz(t, x), Vi/>(i, x), u(t, x)) 
> H(t,x,z(t,x),$(t,x),Vz(t,x),V4>(t,x),v), VveU, a.e. e fl T (3.27) 

and 

H (t, x, z(t, x), xp{t, x),Vz(t, x), Vr/>(t, x), u(t, x)) 

> H(t, x, z(t, x), ■ipit, x),Vz(t, x), V4>{t, x), v) 

+$(A(t, x, v) - A(t, x, u(t, x)), A(t, x, v),Vz(t, x), V^(i, a;), ei), 

Vu e f/, a.e. (t,x) e fir- (3.28) 

Moreover, we can generalize (3.28) to the following: 

i7 (i, x, z(t, x), r/>(i, x), Vz(t, x),Vi>(t, x),u(t, x)) 

> if (f, x, z(t, x), a;), Vz(i, x), Vf/>(i, x), w) 

+$(A(t,x,u) - x, u(t, x)), x, v), Vz(t, x), Vip(t, x), e) , 
Vw e [/, e e S n ~\ a.e. (i, x) e fi T , 

where S" -1 = {x e R" | |x| = 1}. Now, for given (t,x,v), we denote 

A(t,x,v)?e 
^ \A(t,x,v)ie\' 

£ = A(t,x,w)"5 [A(t,x,u) - A(t,x,u(t,x))]Vz(t,x), 
r] = A(t, x, u)~^ [A(i, x, w) - A(t, x, u(f, x))] V^(t, x). 

When e runs over S 71 ' 1 , fi will run over S 71 ' 1 . Then (3.27) and (3.28) become 

H(t, x, z(t, x), i>{t, x), Vz(t, x), V-0(t, x), u(t, x)) 

—H(t, x, x), x), Vz(i, x), Vi/>(i, x), w) 

> sup max (£ T W< T ?7, 0) . 
Ifl=i V 7 



(3.29) 



(3.30) 
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By Lemma 2.5, one can get that 

T T f m±Ah ifn>2 , 

sup £ W V = { 2 

W =1 ^ ifn=l. 



Thus 

l£lM + e T ?? 

p max ^ c; //, u j = 

ImI= 

Combining the above, we obtain (1.9). This completes the proof of Theorem 1.1. □ 



sup max f e w T V, o) = '^1 + ^ (3 31) 

l*il=i V / 2 



We can see that the limit equation (3.10) of homogenizing spike variation equation (3.2) 
may be different for different r. More precisely, (3.10) has essentially three different cases cor- 
responding to r < 2, r = 2 and r > 2, respectively. However, the variational equation (3.18) is 
independent of r > 0. Thus, the final result (Theorem 1.1) can be got by choosing r <E (0,2). 
Such a choice will lead to a simple proof of Theorem 1.1. But, if we did that, we would not know 
whether we can get other conditions from cases of r > 2. This is not satisfied. 

Concerning the method to construct spike variation, we have mentioned that special forms 
of spike variation are needed to get good expressions of the limit equations as (3.10). When 
we introduce (3.1), it is natural to expect that u 4 (-) has no essential effect on the final result. 
The effect oi u 2 (-) is in time scale. One can see that essentially, u 2 {-) works as a spike variation 
as A(-) being independent of u. Difficulties caused by A(-) containing u appear when u 3 (-) is 
introduced. 

Remark 3.1. If we follow the idea of sequential laminates (see Tartar [15]), we can gen- 
eralize (3.10) and (3.17) by constructing more general homogenized equations with their leading 
terms satisfying 

Q s (t,x) = A(t,x) + S(A2(t,x) + A 3 (t,x)-2A(t,x)) 

-S(A 3 (t, x) - A(t, x)) (A 3 (t, x)Y h P (A 3 (t, x))~ h (A 3 (t, x) - A(t, x)) + o(S), 

where 

A(t,x) = A(t,x,u(t,x)), A 2 (t,x) = A(t,x,u 2 (t,x)), A 3 (t,x) = A(t,x,u 3 (t,x)) 

and 

PeYi = {W eS n \W>0,trW = 1} 

n n 

= Q2<x k ZkZk \ai,a 2 ,...,a n e [0, • ■ • ,6. € S , "- 1 ,^a fe = 1}. 
fc=i fe=i 
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Consequently, we can generalize (3.21) — (3.23) with (3.22) being replaced by 

G(t,x) = A 2 (t,x) + A 3 (t,x) - 2A(t,x) 

- (A 3 (t, x) - A(t, x)) (A 3 (t, x))~ h P (A 3 (t, x)y i (A 3 (t, x) - A(t, x)) . (3.32) 



It follows from 



sup £ T Pn = sup !; T xx T r) 

P6./1 |x| = l 



that (3.21), (3.23) and (3.32) still lead to Theorem 1.1. 

4. Problem with State Constraints. In this section, we will consider the cases of state 
constraint. We will only state the results since the proofs are completely similar to those of 
elliptic cases. 

(S6) Let Z be a Banach space with strictly convex dual Z*, F : L 2 (0, T; Hq (Q)) -» Z be 
continuous Frechet differentiable, and E C Z be closed and convex. 

As in Chapter 5 of [10], many state constraints can be stated in the following type: 

F(z{-)) e E. (4.1) 

Let P a d be the set of all pairs (z(-),u(-)) satisfying (1.1) and (4.1). Any (z(-),u(-)) G P a d is 
called an admissible pair. The set U a d = {«(•) € U \ (z(-;u(-)),u(-)) e Pad} is called the set of 
admissible controls. Our optimal control problem with state constraint is 
Problem (SC). Find a control u(-) G U a d such that 

J{u{-))= inf J(u(-)). (4.2) 

To state necessary conditions for optimal controls of Problem (SC), we need to recall the notion 
of finite co-dimensionality (see Ch. 4 of [10], for example). 

Definition 4.1. Let X be a Banach space and Xq be a subspace of X . We say that Xq is 
finite co- dimensional in X if there exist x\, x 2 , ■ ■ ■ , x n G X , such that 

span{A ,xi, • • • ,x n } = the space spanned by {X ,xi 7 ■ ■ ■ ,x n } = X. (4.3) 

A subset S of X is said to be finite co- dimensional in X if for some xq G S, span (S — {xo}) = the 
closed subspace spanned by {x — xq\x G S} is a finite co- dimensional subspace of X and coS = 
the closed convex hull of S — {xq} has a nonempty interior in this subspace. 

Let (£(■), tt(-)) be an optimal pair of Problem (SC). Let Z = Z{-;u{-)) G L 2 (0, T; fl^(fi)) 
be the unique weak solution of the variational equation (3.21) and define the reachable set of 
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variational system (3.21) as 

K={Z(-;u(-)) |u2(-),«s(-)eW}. (4.4) 

Now, let us state the necessary conditions of an optimal control to Problem (SC) as follows: 

Theorem 4.2. Let (S1)-(S6) hold. Let (z(-),u(-)) e V a d be an optimal pair of Problem 
(SC). Let 

F'(z(-))K -E={t- V \te F'{z{-))Tl, r, € E} 

be finite co- dimensional in Z . Then there exists a triple {i>n,i> (•)></?(•)) € RxL 2 (0, T; Hq (fi)) xZ* 
satisfying 

! ^o<0, 

(^oM-))^0, (4-5) 
(^0,^(0) 7^ °) if F'{z{-))* is injective, 
(<p(-),V-F{z(-)))z;z<0, VneE, (4.6) 

ft^t, x) + V- (A(t, x, u(t, x))V$(t, x)) = -ipof%(t, x, z(t, x),u(t, x)) 

-fz{t, x, z(t, x),u(t, x)) ip(t, x) + F'(z(-))*(p, in fi T , 
^(i,x)=0, on[0,T]x9fi, 
</S(T, ar) = 0, in ft, 



(4.7) 



i? (t, x, z(t, x), x), Vz(t, x),Vip(t, x), it(t, x)) — H(t, x, z(t, x),ip(t, x), Vz(t, x),Vil>(t, x),v) 
> ^\A(t,x,v)~%(A(t,x,u(t,x)) - A(t,x,v))Vz(t,x)\ \A(t,x,v)~*(A(t,x,u(t,x)) - A(t,x,v))Vtp(t,x)\ 
+ i (A(t,x,v)~*(A(t,x,u(t,x)) - A(t,x,v))Vz(t,x),A(t,x,v)~^(A(t,x,u(t,x)) - A(t,x,v))Vip(t,x) ), 



where 



\fv e U, a.c. (t,x) e T , 

H(t,x,z,i>,£,r),v) = (ip,f(t,x,z,v))+i) f°(t,x,z,v) - ( A(t,x,u)^,?7), 

(i, x, z, !/),(,i),D)e[0,r]x!lxRxRx R n x R" x U. (4.9) 

Acknowledgement. The author thanks the anonymous referees for their helpful suggestions. 
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